Typeset with jpsj2.cls <ver.l.l> 



Charge Ordering in a-(BEDT-TTF) 2 l3 by synchrotron x-ray diffraction 
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The spatial charge arrangement of a typical quasi-two-dimensional organic conduc- 
tor a-(BEDT-TTF) 2 l3 is revealed by single crystal structure analysis using synchrotron 
radiation. The results show that the horizontal stripe type structure, which was sug- 
gested by mean field theory, is established. We also find the charge disproportion above 
the metal-insulator transition temperature and a significant change in transfer inte- 
grals caused by the phase transition. Our result elucidates the insulating phase of this 
material as a 2kp charge density localization. 
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As an origin of a metal- insulator (M-I) transition, the charge ordering caused by the long range 
Coulomb interaction between electrons is a prominent phenomenon as well as the Peierls transition 
or the Mott Hubbard transition. Organic materials having low-dimensionality provide us a field 
of studying such interesting phenomena. 1 Among them, an organic conductor a-(BEDT-TTF)2i3, 
which undergoes a M-I transition at Tmi = 135 K, 2 attracts considerable attention because of 
its possibility of zero-gap semiconductor under hydrostatic pressures. 3 The insulating phase is a 
nonmagnetic state with a spin gap, 4 and interpreted to be caused by a charge disproportionation 
theoretically 5,6 and experimentally. 7 ' 8 According to these studies, horizontal-charge stripes with 
valences of +le and are expected on BEDT-TTF molecules while no structural evidence has 
been reported. In this letter, we report a precise structure of a-(BEDT-TTF)2l3 as a function of 
temperature based on synchrotron radiation diffraction measurements in order to provide evidence 
of the charge disproportionation. 

a- (BEDT-TTF) 2I3, which has a two-dimensional (2D) conduction band, is a typical organic 
conductor made of I3 ions and BEDT-TTF°' 5+ on average. The crystal structure of this material is 
a sandwich structure of I3 insulating layers and 2D-BEDT-TTF conduction layers having a quarter- 
filled hole band. The space group at room temperature is PI. Together with the metal- insulator 
transition, it shows paramagnetic-nonmagnetic transition. 4 According to several reports, 6, 9,10 this 
transition is not a Peierls transition but a charge ordering caused by a strong Coulomb repulsion. 
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Theoretically, the charge ordering of this material was first discussed by Kino and Fukuyama with 
the Hubbard model, 10 and Seo performed extended Hubbard model calculation 6 that predicted 
several arrangements of charge ordering is stabilized as a function of U /V, where U and V are on- 
site and inter-site Coulomb energy. Experimentally, a so-called "horizontal stripe" type of charge 
ordering is suggested by NMR 7 and Raman scattering 8 measurements. However, no structural 
evidence for this charge ordering has been provided despite of several attempts. 11-13 

Single crystals of a-(BEDT-TTF)2i3 were grown using an electrochemical method. The shape 
of the specimens used for our measurements was a quadratic prism with a typical edge-length of 
0.1mm. Oscillation photographs with 18 keV x-ray for crystal structure analysis were taken with 
an imaging plate (IP) Weissenberg camera installed to beam lines BL-lA and -IB at the Photon 
Factory, KEK, Japan. The sample temperature was controlled by a N2 flow cooler (down to 90 K) 
and a closed-cycle helium refrigerator (down to 20 K). For image data processing of digital IP 
data and structure refinement, RAPID software and CrystalStrucuture software (RIGAKU Co, 
Ltd.) were used. Additionally, we used a four-circle diffractometer having scintillation counter 
installed to BL-4C for symmetry determination based on the Friedel law. The x-ray energy for this 
measurement was 12.4KeV in order to utilize a large value of anomalous scattering factor, and the 
N2 flow cooler was used for controlling the sample temperature. 

Utilizing the Friedel law, the intensity of the Bragg reflection of the Miller index (h, k, I) 
(I(h,k,l)) equals l(h,k,l) for centrosymmetric crystals, we firstly examined the symmetry reduc- 
tion from PI to PI accompanied by the M-I transition, which was suggested by Raman scattering 
measurements. 8 Figure 1 shows the intensity difference SI = for (2, 8, 6) and (2, 8, 6) 

as a function of temperature. The absorption correction was made using the data for higher tem- 
peratures (165K - 180K) at which the crystal is centrosymmetric. The jump in the value of 51 
at Tmi seen in Fig. 1 indicates the breaking of inversion symmetry at the transition temperature. 
Therefore, the space group of insulating phase was found to be PI uniquely. 

In order to make clear the charge ordering structure, we performed crystal structure anal- 
yses at room temperature (RT) and 20K. Lattice parameters at RT were a = 9.1868(13)1, 
b = 10.7926(18)1, c = 17.400(3)1, a = 96.957(4)°, f3 = 97.911(6)°, 7 = 90.795(6)° and 
V = 1695.4(5)1 3 , and at 20K were a = 9.0162(4)A, b = 10.6695(6)A, c = 17.324(1)1, 
a = 96.536(3)°, (3 = 97.758(4)°, 7 = 90.198(4)° and V = 1639.5(1)A 3 . The resolution of the 
data and the number of unique reflections at RT (20K) are 0.48l(0.40l) and 16708 (39714). The 
structure refinement using the full matrix least square method was converged well, showing a good 
reliability factor i? a u = 5.2%(RT) and 5.6%(20K). In the case of an acentric structure, the crys- 
tal has a twin structure of the right-handed coordinated system and the left-handed one, and the 
Flack parameter, which is an indicator of volume fraction between them, was 0.3 at 20 K, meaning 
that our data were taken from a crystal having 70% of the right-handed domain and 30% of the 
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Fig. 1. Intensity difference SI of (2, 8, 6) and (2, 8, 6). The abrupt change at Tmi shows the breaking of the inversion 
symmetry. Solid lines are guides to the eyes. 

left-handed domain. The detailed structural parameters are shown in Table I. 14 The detail of the 
twin structure is discussed later. 

Figure 2 shows the crystal structure of a-(BEDT-TTF)2i3 seen from molecular longitudinal 
direction. At room temperature, there are three crystallographically independent molecules, called 
A, B and C, in a unit cell. Molecule A' is equivalent to A owing to the inversion symmetry. 




Fig. 2. Crystal structure of a-(BEDT-TTF)2l3 seen from the molecular longitudinal direction. Each molecule (A, 
A', B and C) is independent crystallographically below Tmi although A and A' are equivalent above Tmi because 
of the inversion symmetry. The double-side arrows show the labels for the overlap integrals. 

The valence of a BEDT-TTF molecule is estimated by an empirical method based on the bond 
lengths in the molecule. 15 The result of the structure refinement shows the valence conditions of 
A=0.49(3), B=0.57(4) and C=0.41(3) at RT (see upper inset to Fig.3). This charge disproportion 
in the metallic state cannot be caused by the local electric potential distribution originated from 
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Fig. 3. Temperature dependence of the molecular charges. The charge disproportion enlarges significantly at Tml 

the I3 ions, because the valence of the BEDT-TTF molecule that is the closest of ion is 0.4+, 
rather more negatively charged than average. Therefore, the charge disproportion has something 
to do with the transfer integrals among the BEDT-TTF molecules. It is noteworthy that the weak 
charge disproportion was also observed at RT as reported previously. 8 ' 16, 17 

In the low-temperature insulating phase, we have succeeded in a direct determination of the 
charge ordering structure using the same method as above. The space group is PI and there are 
four independent molecules in the unit cell, or molecule A' is independent on A. We observed 
molecular deformations, which reflect charge disproportion, at 20K. Estimated each charge value 
is A=0.82(9), A'=0.29(9), B=0.73(9) and C=0.26(9). The charge structure is a horizontal stripe 
type, in which molecules A and B are hole rich and A' and C are hole poor sites as shown in the 
lower inset to Fig. 3. The obtained arrangement corresponds to the theoretical predicted one by 
Seo. 6 

To investigate the electronic and structural change associated with the charge ordering tran- 
sition, we performed the structure analyses at intermediate temperatures, 210K, 175K, 147K and 
110K. The reliability factors of structure refinements are small enough (ca. 5%) for making a 
precise discussion of electronic state. Figure3 shows the temperature dependence of the molecular 
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charges. The relation of the molecular charge (B>A>C) observed at RT conserves in the metallic 
state, while the charge distribution changes gradually. Based on the tight binding approximation 
and a molecular orbital calculation with the extended Huckel method, 18 we calculated the overlap 
integrals \S\'s between the neighboring molecules, which are proportional to the transfer integrals 
approximately. Figure 4 shows the temperature dependence of the overlap integrals. Labels in 
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Fig. 4. Temperature dependence of the overlap integrals. Because of symmetry reduction from PI to PI, there are 
twelve integrals below Tmi although seven integrals in metallic phase. All labels correspond to those in Fig. 2. 

the figure correspond to those of Fig. 2, in which the \S\'s having labels with and without an 
apostrophe are equivalent in metallic phase because of the inversion symmetry. The number of the 
overlap integrals increases from seven in the metallic phase to twelve in the insulating phase due 
to the reduction of symmetry. 

It is noteworthy that the overlap integrals are able to be classified into two groups; bl and b2 
form a main transfer path regarded as a one dimensional "zigzag" chain in the unit cell and the 
other paths form inter-chain interaction as shown in Fig.5(upper). The values of {Si's for bl and 
b2 in the chain increase with decreasing temperature down to Tmi- 

Below Tmi, overlap integral (S^bl which corresponds to the hopping from molecule A' to 
molecule C become fairly small and the values of \S\h2 and |<S'|bl / approach each other. The 
modulated transfer integrals in a zigzag chain are schematically shown in Fig. 5(lower). The 
structure is very similar to the 2&p charge ordering, where kF is the Fermi wave number. This 2Uf 
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Fig. 5. (left) Schematic view of a zigzag chain made of large overlap integrals |5*|bl and |S|b2. Thickness of the line 
shows the magnitude of the overlap integrals, (right) Expanded picture from zigzag chain to a one dimensional 
chain. In the low temperature phase, there are two possibilities of the valence arrangements to take. 

structure causes both charge gap and spin gap, and makes the system a nonmagnetic spin singlet 
state. From these results, we can recognize that the electronic (band) structure exhibit an abrupt 
change accompanied by M-I transition. 

The main structural change creating this modulation of transfer integrals is not translational 
shifts of molecule but changes in the dihedral angles. The maximum degree of change in the 
angles between RT and 20K was about 3°, and the change in the angle make the spatial overlap 
of the wavefunctions. This structural charge ordering pattern would be described by the ordinary 
low-dimensional electronic state and is in striking contrast to Wigner crystallization of the charge 
ordering in quarter filled conductor, (DI-DCNQI^Ag system. 19,20 

Finally, we mention the twin structure. When the crystal undergoes the phase transition from 
PI to PI, there are two kinds of structural deformation. The situation for this particular case is 
depicted in Fig. 5. One side of the twin is the (A and B)-rich type and the other is (A' and B)-rich 
type. The left-side panel for the low-temperature phase in this figure shows the former. The overlap 
integrals for bl' and b2' are larger than those for bl and b2, respectively, in the former, and the 
other way around in the latter. The differences between them are only the phase of charge density, 
and the electronic states of them are completely equivalent. 

We succeeded in the determination of the charge arrangement in the insulating phase of a- 
(BEDT-TTF)2l3 by means of synchrotron radiation x-ray crystal structure analysis. Based on the 
molecular deformation, we found horizontal type of charge ordering directly. In the metallic phase, 
charge disproportion in the conduction plane was also observed. The result of a quantum chemical 
calculation based on the obtained structure parameters shows a drastic change in transfer integrals 
accompanied by M-I transition. The insulator state is described with 2kp charge density wave 
within a zigzag path connected by the large transfer integrals. 
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Table I. Refined fractional coordinates and isotropic displacement parameters at 20K. 
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